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I Introduction 



C I The tunneling of particle through potential barrier is an essentially quantum phe- 
^ ' nomenon. This process involves propagation of a particle through a classically inac- 
O^' cessible region. The complete information about tunneling of particle from solution 
of Schrodinger equation with appropriate boundary conditions can be obtained in 
^ \ all regions. But in practice, the exact solution of Schrodinger equation can be found 
\ for some simplest forms of potentials and it is difficult enough to obtain the exact 
solution for arbitrary potential form. For this reason the approximation methods are 
used for finding solutions for potentials of specific form. But exact solutions, which 
were found, are of a great importance, because they allow to analyze the tunneling 
process in general. 

It is enough difficult to obtain solutions for multi-dimension potential forms. 
Therefore, in present work only one-dimensional case is considered, for which exact 
solutions are obtained for some simple forms of potential having two wells separated 
by barrier. Every specific form of potential in comparison with another ones requires 
use of specific approach for solving the problem, allowing some features of tunnel- 
ing to be more pronounced, while the others will be more unnoticeable. Having 
analyzed the use of boundary conditions for solving the problem, we proposed to 
divide various shapes of double-well potentials into two classes: squared potentials 
and potentials having rounded off forms. Squared double-well potentials have exact 
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analytic solutions (which can be expressed through elementary functions). Rounded 
off double- well potentials have solutions expressed through special functions (if these 
solutions exist). In present work after qualitative analysis of energy levels of various 
forms of double-well potential the problem of squared potential and the problem of 
rounded off potential (the Morse's potential) are considered separately in next two 
sections. 

The period of particle oscillations between two wells is one of the most important 
parameters which characterize the process of tunneling. We can obtain the value of 
a period on the basis of energy levels of system. For this at the beginning we will 
confine ourselves to class of systems, for which the distances between energy levels 
have the exactly determined largest common devizor A and the following condition 
is fulfilled: 



En = Eo + A* In, (1.1) 

where In & 0,N. In general case in region of discrete energy spectrum the states 
of these systems are described by wave packet as follows [T] 

m =Y.9n^n{x)e-''^''--'''^"'' = Y.gnVn{x)e-''''""'' (1.2) 

n n 

where ({>n{x) is orthonormal wave functions of stationary states of system satis- 
fying the equation H '^n{x) = En'-Pnix)'-, H is Hamiltonian of system; Ifl'nP = 1 

n 

and insignificant factor e"*^"*/^ is omitted, which is common for all terms of sum 
J2- Let's select the moment t = as the origin of time reference. 

n 

Let the function ijj{t) be determined on region [— ^] and satisfy the Dirichlet's 
conditions [2]: a) it can divide this region into finite number of regions, in which the 
function ip{t) will be continuous, monotonic and bounded; b) if to is the discontinuity 
point of function ip{t), then ■^(to + 0) and '?/'(^o ~ 0) exist. Then expression (1.2) is 
expansion in t of function ip{t) into Fourier series which is converging in all points of 
region [— Then the function ■ip{t) is periodic in time and period of oscillations 
(time of Poincare's cycle) is given by 

2TTn , , 

^ = ^ (1-3) 

The expression (1.3) determines period of oscillations of wave packet, if energy 
levels in region of discrete spectrum have exactly determined common largest divizor 
A and condition (1.1) is fulfilled. 

In general case for systems, for which the distances between energy levels in 
region of discrete spectrum don't have the exactly determined divizor A, one can 
select the 'quasi-cycles' with given degree of accuracy, for which the state of system 
approaches the maximum degree the initial state after a 'quasi-cycle' time |3j . States 
of such systems are localized in a confined volume of space and the time of Poincare's 
cycle (which includes the required number of 'quasi-cycles') can be determined with 
given degree of accuracy. To find information about some parameters for quantum 
systems evolving with time in region of discrete spectrum see |Sl IH E] • 

Since period of particle oscillations between wells is obtained on the basis of 
energy eigenvalues, much attention is paid to the problem of solving of eigenvalue 
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equations. For some forms of potential the transmission coefficient through the 
barrier is found. This parameter can be obtained using two approaches: in region of 
continuous spectrum for particle which is incident upon the barrier, with asymptotic 
velocity hk/m, and, when condition of WKB approximation is fulfilled, in region 
of discrete spectrum for particle which initially is located in one well and then is 
tunneling to another well. In case of double-well squared potential the comparison 
of these two approaches for calculating the transmission coefficient through barrier 
is performed. For symmetric double-well potentials the dependence of transmission 
coefficient (which is found using one of approaches) on period of particle oscillations 
between wells is analyzed. 

II Analysis of possibility of particle tunneling through 
barrier 

The qualitative estimation — will the particle oscillate between wells or not — can 
be obtained on the basis of solution analysis of stationary Schrodinger equation 
which are found for every well. 

Divide various forms of double-well potentials into two classes: squared potentials 
(see Fig. Hj) and potentials, which have rounded off forms (see FigE)). To obtain the 
energy levels, the Schrodinger equation is solved separately in every region (3 regions 
for squared potential, see Fig.P 2 regions for potential having rounded off form, see 



As a result of this solution, the wave functions ipi{x) are found in every region, 
and these functions must be continuous and bounded (we consider discrete spectrum 
E < Uo) all over the region of its determination (here, i is a number of region). If 
wave function is expressed through special or elementary functions, then it will 
be bounded everywhere in region —d < x < b (with exception of some cases of 
hypergeometric functions, which must be considered separately, see 0111), and in 
points X = —d and x = h the boundary condition is given by 



This condition determines eigenvalues for energy levels as for discrete levels. 
Condition of continuity of wave function and its derivative in region —d<x<b 
requires the equality of solutions ip{x) and dip{x)/dx for adjacent regions in points 
of boundary between these regions (in these points the discontinuity of derivative is 
possible) . 

At the beginning we consider squared potential (see Fig.lH). For particle located 
in left well (region 1) we analyze the possible cases of its propagation to right well 
in result of tunneling through barrier. Among them we select the following cases: 

1) In regions 1, 2 and 3 there are the equal levels Eq (see Fig.Q), which are found 
from Schrodinger equation for every region separately. Then particle can propagate 
through barrier along level Eq. (In this case the transitions between levels are not 
required for tunneling). We obtain the eigenvalue of energy level Eq from system 
(2.1) and the following system: 



Figiai. 
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<^3(a) = Ma), M(^) = ^2(a)- 

2) In regions 1 and 3 there are the equal levels Ei, but such level is absent in 
region 2. In regions 2 and 3 there are the equal levels E[, but such level is absent in 
region 1. Then particle initially located on level Ei in left well, can not propagate 
to right well along this level. But in region 3 wave functions corresponding to levels 
El and E[ are not equal to zero. Therefore, in this region the matrix element of 
transition from level Ei to level E[ (and on the contrary) is not equal to zero. 
Therefore, particle which initially located on level Ei in left well can propagate to 
right well with transition from level Ei to level E[. The transition Ei E[ takes 
place in region of barrier 3. Eigenvalue of level Ei can be obtained from the system: 

(^l(-c) = ^3(-c), <^i(-c) = M-C), ,^ „^ 

^iM)=o, ^3(a) = o. ^^-^^ 

In this fashion, one can find the eigenvalue for level E'^ from the following system: 

<^2(a) = <^3(a), V52(fe)=0, ^2 4) 

'^2(0) = V'sW- ^3(-c) =0. 

3) In region 1 the particle is located on level E2 and this level is absent in regions 
2 and 3. Also in regions 1, 2 and 3 there are the equal levels E'^. On this case the 
particle can not propagate along level E2 to right well and to region 3 of barrier 
(there is a full reflection of particle along level E2). But the particle can make 
transition on level E2I in region 1 and then it can propagate to regions 2 and 3 along 
this level. One can find the eigenvalue of level E^ from the following system: 

V9i(-d) = 0, ^i(-c)=0. (2.5) 

The eigenvalue of level E2I satisfies system (2.1) and system (2.2). 

More concretely the first case will be considered in one of the next sections. 
Note, that at symmetry of potential {d = b, c = a, Wq = 0) for every level of left 
well one can find the appropriate level in right well and on the contrary. Therefore, 
for cases considered above only the first case is possible for this potential. 

Now we consider potential which has rounded off form (see Fig. 12)). Let the 
particle be localized in the left well. Among possible cases, in which the particle 
can propagate through barrier to right well, we select the following cases: 

1) In regions 1 and 2 there are the equal levels Eq. Then particle can tunnel 
from left well to right well along level Eq. The eigenvalue for this level can be found 
from system (2.1) and the following system: 

MO)=MO), MO) = MO)- (2.6) 

2) In left well the particle is located on level Ei and this level is absent in region 
2. But in regions 1 and 2 there are the equal levels E[. Then particle can not 
tunnel from the left well into the right one along level Ei. But at first it can make 
transition from level Ei to level E[ in region 1 and then it can propagate to region 
2 along level E[ . System for finding eigenvalue of level Ei is given by 
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^,{-d) = 0, </^i(0) = 0. 



(2.7) 



Eigenvalue for level E[ can be obtained from system (2.1) and system (2.6). 

In one of the next sections the problem with potential of such form will be 
considered more concretely. As an example, the double-well Morse's potential is 
selected. Also note, that for potential of rounded off form (as in case of the squared 
potential) only the first case is possible when the potential is symmetric. 

Ill Dependence of distance between two closely 
located levels on transmission coefficient 
through barrier in quasi-classic symmetric po- 
tential 

Consider potential U{x), which has two symmetric wells separated by barrier (see 
Fig. El). If barrier is not penetrable, then there are energy levels Eq corresponding 
to oscillations of particle only in one well. The possibility of particle transitions 
between wells leads to splitting of every level Eq into two closely located levels Ei 
and E2. 

We consider the case, when potential U{x) is quasi-classic. Then the splitting 
value AE can be obtained through wave function (po{x), which determined with 
accuracy to first order terms in ^, as follows jH] 

AE = E2-E^ = ^<^o(0)(^[,(0) = — * exp (-\ J \p\dx) , (3.1) 

—a 

where 



^o(O) = 




and Vq = \j2{Uo — E)/m, p is impulse of system, w = 27r/T is frequency of 
classic periodic oscillations; a is turning point corresponding to level Eq (see Fig.lHI). 

Transmission coefficient D through barrier in WKB approximation is determined 
by i 



D = const * exp 




—a 



where proportionality factor const is determined by accuracy of WKB approx- 
imation and is equal to 1 with accuracy of first order terms in h (see jH]). Taking 
the expressions (3.1) and (3.3) into consideration, we can write 

AE = const * — (3.4) 

TT 
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Here, transmission coefficient D is determined by expression (3.3) for discrete 
energy spectrum. In accordance with conditions of using WKB approximation, the 
expression (3.4) determining the dependence of level splitting AE on transmission 
coefficient D is used only for small values D. 

Now we consider some cases, for which there are the exact analytical solutions. 



IV Double-well infinite squared potential 

Consider system, potential of which consists from two squared wells separated by 
squared barrier of finite height (see Fig. Q). This potential is given by 

cxD, for X < —d, X > b] 

0, ioT — d < X < — c, (region I); 

f/o, for — c < X < a, (region III); 

—Wo, for a < X < b, (region II). 

In case of discrete energy spectrum in region Uq > E > we find the solution of 
stationary Schrodinger equation in form: 

{Oi sm{k{x + d)), 
as sin(fc3(x - 6)), 
Here, the following coefficients are used: 

k = lV2mE, 

= y2m{E + Wo), (4.3) 
X = \pm{Uo-E). 

Let's consider the case of particle oscillations between wells along one energy 
level (without transitions between energy levels). Unknown coefficients and these 
energy levels can be found from continuity conditions of wave function in boundary 
points x = — c, X = a and from the following normalization condition: 

+00 

J \ip{x)\'^dx = 1. (4.4) 
—00 

In result, we obtain unknown coefficients: 




for — d < X < — c, 

for - c < X < a, (4.2) 
for a < x < 6. 



ai 



d — c a + c 



sin^^k^d — c)) cos'^{k{d — c)) 



X 



[ 2 2 

/ sm{k{d-c))-k/xcos{k{d-c)) \ 2 ^ ^_2x(a+c) 
\sin(A;3(6 — a)) + ks/xcos{k3{b — a)) J 

* sin(2A;3(a — b)) 



+ 1 



b — a 1 

\ * 

2 4A;3 



(e2x(a+c) _ 1) + 



\ sin(2A;(c^ - c)) + (M^jd - c)) + k/xcos{k{d - c))f ^ 
Ak 8x 

-1/2 



{sai{k{d — c)) — k/xcos{k{d — c))Y 
8^^ 



(4.5) 
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02 = cti * ^e^'^{sin{k{d — c)) + k/xcos{k{d — c))), 
62 = cti * |e~^'^(sin(/c((i — c)) — k/xcos{k{d — c))), 



03 = oi * e 



^(a+c) fe/ X cos{k{d - c)) - sm(A;(t/ - c)) 
ks/xcos{ks{b - a)) + sin(A;3(6 - a)) ' 
Eigenvalue equation for this potential are given by 



(4.6) 



E = h^k^/2m, 

sin(fc3(6 — a)) — k^/xcos{k^{h 



a)) sm{k{d — c)) — k/xcos{k{d — c)) 



^ g2x(a+c) _ 

sin(/c3(6 — a)) + k^/xcos{k2,{h — a)) sva.{k{d — c)) + k/xcos{k{d — c)) 

(4.7) 

Now we consider symmetric case {d = b,c = a,Wo = 0) jHJin]. Wave function 
became symmetric or antisymmetric: 



(4.8) 



ai sin(/c„(x + 6)), 
^n{x) = <; 62((-l)"e>^" + e-'^^), 

ai(— 1)" sin(A;„(6 — x)), 

where n E 0,N. Unknown coefficients ai and 62, obtained from expressions (4.5) 
and (4.6), are given by 



for — b < X < —a; 
for — a < X < a; 
for a < X < b. 



(sin(fc(6 - g)) - fc/x cos(fc(& - a)))^ ^ , ^ . 
^ - « + ^^55^^ * I (-1) 2a + 

sin(/i;(6 — a)) — k/xcos{k{b — a)) 
The wave vector k is transformed to form: 




(4.9) 



k 



1 



* < — arcsm 



'1 + 



'_l'jn _ g2xa^ 2 



+ nn 



(4.10) 



P V(-l)'' + e2x«. 

Consider the particle propagating from the left to right in the potential (4.1) 
with asymptotic velocity hk/m and energy E < Uq. For it the wave function can be 
written as follows 



-1/2 ^ l^^ikx _ ikx-^^ ^Qj^ ^ ^ _g. 

'1/2 * (5e>^^ + E'e-^""), for - c < x < a; 
"-■^/^ * Ce^^^, for X > a. 



(4.11) 



{(/iA;/m) 
(/;,A;/m) 
(/i/c/m) 

Coefficients A, 5, 5' and C are obtained from continuity conditions for (p{x) and 
dip{x)/dx at points x = — c and x = a. The transmission coefficient D and reflection 
coefficient R calculated as the ratio of the flux of incident wave in region III to the 
flux of transmitted wave in region I or reflected wave in region III, are given by 



D 



R 



Akksx' 



{kk^ + x'^y sinh2(x(a + c)) + x^{k - k^Ych^{x{ci + c)) + Akk^x^ 

{kks + x^f sinh2(x(a + c)) ^x^ik- k^fch^{x{a + c)) 
{kk^ + x^Y sinh2(x(a + c)) + x\k ' k3)W{x{a + c)) + 4fcfc3x2 ' 



(4.12) 



We find the values D and R for transmission of particle through barrier in con- 
tinuous energy spectrum. But comparison of transmission coefficient D with its 
small values determined by expression (4.11) for symmetric case with transmission 
coefficient D determined by expression (3.3), which is obtained in WKB approx- 
imation for discrete energy spectrum, show, that both approaches give identical 
formulations with accuracy to normalized constant: D = const * exp(— 2xa), where 
const is determined by accuracy of WKB approximation and is equal to 1 in terms 
of first order in [SI In this sense, we will formally consider the expression 
(4.11) as the determination of transmission and reflection coefficients for discrete 
energy spectrum. The values k and X; used in expression (4.11), can be obtained 
from equation (4.7) or (4.10). 

We analyze the periodicity of particle oscillations between wells in symmetric 
potential. For this we consider equation (4.10) which can be written as 

fi{kn) = kn{b-a), 

f2{kn) = -arcsin| ^ .,2 ./^^n _ g2xa. 2 } + (^'IS) 




-1)" + e2x«/ 

The graphic analysis of exact solutions of system (4.13) gives a number of values 
k^ (see Fig. IH). Changing the second equation of system (4.13) to its linear relation 
f2{kn), which can be obtained by linear approximation, we flnd the following values 
of kn- 

kn = ko{2n + 1), 

Er, = Eo{2n + 1)2 = Eo + AEq * n(n + 1) = + A * (4.14) 
A = 4 * Eq. 

where n, /„ G 0,A^. From expressions (4.14) we obtain period T of particle 
oscillations between wells: 

T = ^ = ^ = -*^^^^ = |^*(2n + l)2. (4.15) 



A 2Eq 2 E„ hk, 



For values kn the accuracy 7r/(2(6 — a)) is used. The series of solutions /c„ is 
bounded by maximum value kj^, where N (the number of energy levels in region 
E < Uq) satisfles the following condition: 



^ V2m[/o - hkp 
2hko 

From expression (4.15) we obtain: 

i; = (2„+l)^- = A*^^. (4.16) 

Using expression (4.12) for symmetric case and expression (4.16), we flnd the 
dependence of transmission coefficient D through barrier on oscillation period T 
and on value A: 



8 



2, /2mt/o ^2^^^ 



sinh^(2aW— ^-(2n + l) , ^ 

DniT) = \l + 1^ . (4.17) 

1 4(2n+l)^^(l-(2. + l)^^)/ 



,2, 2mUo , (2n+l)^ 
sinh^(2aW — ^ - A * ^ ^ 



r V 4f/o 

Also note, that in both hmiting cases D and D —* oo the oscillation period 
T approaches a value: 

T=l!^^^^ (4.19) 
Trn. 

V Double-well Morse's potential 

Consider the system, potential of which is shown in Fig. |21 and is given by 

TK \ - j A(e-2"(^+'=) - 2e-"(^+'=)), for - < x < 0, , . 

~ \ 5(e2/^("-'') - 2e^("-")), for < X < 6. ^ 

For this potential we find solutions of stationary Schrodinger equation. Perform- 
ing the following changes: 



= ^^^e-°(^'+^), for -d<x<0, 
ah 

jjh 



(5.2) 



and introducing the following notations: 



V2mA ( 1 

sa = 1 , nA = Ua + tt 

ah \ 1 



V- 


-2mE 


ah, 


V- 


-2mE 



we transform the stationary Schrodinger equation to the following form: 



(5.3) 



(iV '^d'^ / I n + s + 1/2 



+ 737 + -7 + - 7^ 



= 0. (5.4) 

A,B 



Here, the indexes A or i? are used by values (p, ^, n and s with dependence 
on that, the equation (5.4) is solved in region of left or right well, respectively. For 
obtaining the general solution, let's omit these indexes. Using the following changes: 

c = l + 2s, ^^-^^ 
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the equation (5.4) can be transformed to confluent hypergeometric equation as 
follows 



The particular solutions of this equation can be written by use of confluent 
hypergeometric function F(— n, c;^) jElIZ!: 

1/2(0 = 0-^F(-n-c+l,2-c;0, 7^ 
1/3(0 = e^F{c + n,c--0, ^ ' 

2/4(0 = 0-^e«F(l + n,2-c;-0. 

Consider the case c ^ Z (i.e. 2s ^ Z). Then both the particular solutions and 
1/4 are linear dependent on yi and 1/2, because they are transformed to yi and 1/2 by 
Kummer's transformation [HI Ej : 

F(a,c;0 = e«F(c-a,c;-0. (5.8) 

Solutions ?/i and ?/2 are linear independent between themselves. It can see this 
from its behavior close hj C, 0. Then the general solution of equation (5.6) can 
be written as follows 

1/(0 = CiF(-n, c; + C20~'=i^(-^ - c + 1, 2 - c; (5.9) 

where ci and C2 are arbitrary constants. In initial variables the general solution 
has form: 



^(0 = e-^c^Fi-n, 1 + 2s- + C2r'F{-n - 2s, 1 - 2s; 0] (5.10) 

Expression (5.10) with relations (5.2) and (5.3) is general solution of stationary 
Schrodinger equation at 2s ^ Z and E < for regions of left and right wells. 
Coefficients cia, cib, C2a and C2b (where indexes A 01 B for coefficients ci and C2 
indicate at left or right well, respectively) are obtained from boundary conditions 
and normalization condition. Imposed boundary conditions determine the energy 
spectrum as a discrete ones and using them it can find the energy eigenvalues of 
system. 

Discrete energy spectrum determines the requirement that wave function will be 
finited all over the region of its determination. Finity of wave function in regions 
—d < a; < and < s < 6 depends on constraining condition of particular solutions, 
by use of which the general solution (5.10) can be represented. Constraining of 
particular solutions at finite values x depends on constraining of functions F{—n, 1 + 
2s; and F{—n — 2s, 1 — 2s; 0- But these functions are confiuent hypergeometric 
and are represented by converging series at finite ^[HEj- Atn = OornGA^ the 
first function has form of polynomial, and therefore, it is bounded (because values 
^ are finited). In this fashion, at n + 2s = or n + 2s G the second function is 
bounded. Therefore, both the particular solutions are bounded all over the region 
—d < X < b (factor at ^ — > is finited because of constraining condition at 
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X = 0). Therefore, the general solution (5.10) is also bounded all over the region 
— oo < X < oo. 

Since wave function satisfies condition of constraining all over the region of def- 
inition X, it can be normalized by 



u 

J \ip{x)\'^dx ^ 1. 



(5.ii; 



Now we find the energy eigenvalues of system. Under consideration of possible 
solutions we select two cases: 

1) The particle oscillate between wells along one energy level (without transitions 
between energy levels). On this case the energy levels, which are obtained from 
eigenvalue equation in region of left well, must correspond to energy levels, which 
are obtained from eigenvalue equation in region of right well. On this case the 
following boundary conditions are required: 



(5.12) 



(pA{-d) = 0, 
^B{b) = 0, 

where <fA{x) and ^Pb{x) are the general solutions (5.10) of stationary Schrodinger 
equation in regions of left and right wells, respectively. Solution of equation system 
(5.12) gives the eigenvalue equation of form: 



'Pi 



BO 



( CiA 
*< * 

lC2A 



— * ^BlFi-nB, 1 + 2sb; ^bo) + ^m^Fi-riB - 2sb, 1 - 2sb; i 

C2B 

1 



BO 



^tAoFi-UA, 1 + 2sa; Uo) + sa * Go 'F{-nA, 1 + 2sa; Uo)+ 
+ 



^,^ dF{-nA,l + 2sA;Uo) ] , r l^-.AW „ 9. 1 9. ^ 



AO I 



SA-i-p( o 1 o t ^ , ^-SA ^F{-nA-2sA,l-2sA]i 
-SA * Uo F{-nA - 2sA, 1 - 2s A] Uo) + Uo ^7 

= {— * \-\imF{^nB. 1 + 2sb\ ^bo) + sb * Cbo'^^I-^b, 1 + 2^^; ^bo) + 

I C2B L Z 

,^ dF{-nB,l + 2sB\iBo) ] \ ^c-SBT^i ^ 9« 1 9. ^ 

+?B0 ^7 + -l^^BO ^\~^B- ^SB.i-- ^Sb.Kbo)- 

Oi;B J L Z 

.-SB-iTpf o 1 o t N I ^-SB dF{-nB - 2sb, 1 - 2sb; j 

-SB * ^Bo F{-nB - 2sb, 1 - 2s b; ^bo) + Cbo ^7 

Cc,b 

* [— * e/oFi-riA, 1 + 2sa; Uo) + Go^Fi-riA - 2s a, 1 - 2s a; i 

IC2A 



Uo 



(5.13) 



where 



CiA 

C2A 
ClB 

C2B 



ad\-^'^ F{-nA-2sA,l-2sA;Uoe'"') 
) * Fi-nA,l + 2sA;Uoe-') ' 
fc pb\-^'B F{-nB - 2sb, 1 - 2sb; ^bo^'^^) 
^ * F{-nB,l + 2sB;iBoeP') ' 



(5.14) 
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2V2mA 
? AO = z — * e 



—ac 



Transform expressions (5.3) to form: 




t(^a + «a + - 



/a 1\ 



To find the eigenvalues, we need to substitute the expressions (5.14) and (5.16) 
into equation (5.13) and to resolve it relatively the value sa-, which unequivocally 
defines eigenvalue E. The solution of equation (5.13) is performed by use of nu- 
merical methods and determines the energy eigenvalues corresponding to particle 
oscillations between wells. 

2) Now consider another case: the particle oscillates in one well (for example, 
in left one). On this case its full refiection take place from the middle of barrier 
(here, the transition of particle to another energy level is possible, which exist in 
both regions, with further tunneling of particle along it). Note, that on this case 
the refiection of particle from barrier is principally possible along energy levels of 
range Q > E > Uq (this case of refiection is impossible for symmetric potential). We 
determine the condition of particle refiection from barrier as follows 

(^^(0) = 0. (5.17) 
Using this condition and also the following boundary condition 

^A{-d) = 0. (5.18) 

we obtain equation, from which it can find the energy eigenvalues for particle 
oscillations in left well: 



^2s^ad ^^(_^^ _ 2s^, 1 _ 2sA] Uo) * F{-nA, 1 + 2s a; ^Aoe"'') = 

= F{-nA - 2s A, 1 - 2s A] ^Aoe"'') * F{-nA, 1 + 2sa; Uo) (5.19) 

Resolving this equation relatively unknown values s^, it can find the energy 
eigenvalues of system. In this fashion, it can obtain equation which determines the 
energy eigenvalues for particle oscillations in right well: 



e^sBfib ^p^_^^ _ 2sj,, 1 - 2sb; iso) * F{-nB, 1 + 2sb; ^Boef"') = 

- F{-nB-2sB:l-2sB;iBoe''')*F{-nBA + '2sB;iBo) (5.20) 

Equations (5.19) and (5.20) include the case when some energy levels of left well 
are equal to some energy levels of right well. On this case the following condition is 
fulfiUed: 
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V?a(0) = ipB{Q) = 0, 



^AiO) = ^b(O) = 0, 

which correspond to particle oscillations between wells. Therefore, it need to 
except such energy levels from analysis of particle behavior in one well. 

Consider symmetric potential. On this case every energy level obtained by so- 
lution of eigenvalue equation in region of left well, is equal to corresponding energy 
level which is obtained by solution of eigenvalue equation in region of right well. 
Therefore, the particle located on any energy level of region E < will be oscil- 
lated between wells. Wave function becomes symmetric or antisymmetric. System 
of equations for finding energy eigenvalues is given by 

if wave function is symmetric (even states), and 

^{-d)=0,<f'{0)=0, 
if wave function is antisymmetric (odd states). 

VI Symmetric potential of form + B'^/x'^ 

To analyze the particle behavior in symmetric potential with enough high barrier, 
it can use potential of the following form: 

f/(x) = ^(x2 + — j (6.1) 
where x g] — oo; +oo[. This potential is shown on Fig. El Use new parameters: 

2mE 

F = (6.2) 

K = TT- * B\ 

Then stationary Schrodinger equation is transformed to form: 

^+fG + Fx2 + :^)*¥. = (6.3) 
ax^ \ x'^ ) 

Let's find solutions of this equation. Perform the following change of variables: 

^ = ax^ , 

mw I — — 

a = — = v^, (6.4) 
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In result, the equation (6.3) is transformed to standard Whittaker's form [HIIZ]: 



+ 



1 G 1 K\l_ 

■4^ 4v^e ^ Vie ^TJe" 



w = (6.5) 



Using the following parameters and performing the following change of variables: 



k 



G 



,2 



4v^' 
1 K 



^ - J6"T' (6-6) 
a = - - /c + /i, 
c = 1 + 2/i, 



^(0 = r^/v/2^(0, (6.7) 

we transform the equation (6.5) to confluent hypergeometric equation of form: 

Particular solutions of this equation can be represented by confluent hypergeo- 
metric function F(a,c;^) as follows 

1/1 (0 = F(a,c;0, 

?/2(0 = e-''F{a-c+l,2-c;0, q^ 

7/3(0 = e«F(c-a,c;-0, ^^'^^ 

1/4(0 = 0-'=e«F(l-a,2-c;-0. 

Let's consider the case c ^ Z. In accordance with Kummer's transformation 
[HI Ej, the solutions 2/3 and 7/4 can be written through yi and 1/2. Therefore, the 
solutions 1/3 and 7/4 are linear depended on yi and ?/2- Write flrst two solutions yi 
and y2 in initial variables: 

(x) = ai/2+Ma,-i+2/.g-axV2^ (^^^ I ^ 2/i; ax^) , 
<^2(x) = ai/2+Ma;-ie-°"'/2F (a - 2/i, 1 - 2/i; ax^) . ^ ^ ^ 

Condition of wave function flnity all over the range x requires, that the following 
conditions are fulfilled: 

for V9i(x) : a G 0, -A^; 2/i ^ -A^, ffi 111 

for ^2(a;) : -a + 2/i G 0,A^;2/i ^iV. ^ ' 

As a result of these conditions, the series, which represent confiuent hypergeo- 
metric function for solutions (pi and <^2, transform to polynomial and the spectrum 
becomes discrete. Analysis of expressions (6.11) shows, that solutions ipi and ip2 
can not be used at same time. But both solutions have equal energy eigenvalues E 
which can be written as 



£j = 2A.(i+„±l/l + l!^fl=) (6.12) 
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where n G 0, A^. The general sohition of wave function can be represented 
through ipi or ip2 (which are hnear independcd between themselves). Both energy 
eigenvalues and E~ used in expression (6.12) can be considered separately as 
defining two independed waves with oscillation periods and T~ , respectively. To 
obtain period for every wave, we write the equation (6.12) as follows 



^E^ + A*n^ 2hw *n + 2hw * ( - ± -\/l + 



From this we obtain 



1 1 / Am^w^ 



Et = 2r™.(-±-^l + — 
A = 2hw, (6-13) 



2TTh TT 



A w 

It can see from expressions (6.13), that two considered waves have equal periods. 
In accordance with the following relation 

it can see that general period of oscillations between wells is equal to or T": 

T^T^^-. (6.14) 
w 

Now we find dependence of transmission coefficient D through barrier on oscilla- 
tion period T between wells. Consider the case of small values D, when expression 
(3.4) can be used. From expression (6.13) we find the distance between two closely 
located levels by 



= - E- =hw^l + ^^B^ (6.15) 

Using this relation and also expression (3.4), we obtain the transmission coeffi- 
cient as follows 



4m^ 2^ 2 A I m'^B^ 
D r. n (l + -^B)=n [1 + -^^ 

This expression establishes the dependencies of transmission coefficient D on 
oscillation period T between wells and on largest common divizor A which is deter- 
mined by expression (1.1). 
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VII Double-well symmetric parabolic potential 



Consider the system, potential of which can be written as 

. f \m'u?(x ^ of , for X < 0; . . 

^ ^ 1 ^mw [x — ay, for X > 0. ^ ^ 

Let's assume, that potential U{x) satisfies conditions of using WKB approxi- 
mation. In result of tunneling through barrier the displacements of energy levels, 
which takes place because of level splitting, from their positions without tunneling 
are determined by expression (3.1). For potential (7.1) the energy eigenvalues are 
given by 

E- = hw{n + I) - AEn, 

E+ = hwin + i) + AEn, /7 n) 

where E~ and E^ are eigenvalues, which occur because of splitting and corre- 
spond to symmetric and antisymmetric wave function, respectively. Lp^^ is solution 
of stationary Schrodinger equation in region x > without splitting and has form: 

<^i°)(x) = A„e-"'(^-'^)'/2^„(a(x - a)) (7.3) 



where a = ymw/h, Hn{^) is Hermitian polynomial. Using (7.3), we find the 
displacement value: 

fl^ 2 2 

AE„ = —Ale-'' " HJaaMa^aHJaa) - 2nHn-i(aa)) (7.4) 
m 

Coefficient An can be obtained from normalization condition of (f^\x). Function 
ip!^^{x) is assumed as normalized one, that integral of Iv^^^-'P in region of right well 
is equal to 1. Then it can write 



oo 

2/^_„\2 



Al = \ I e"" iy^(«(x-a))rfx| 

^(^(1 + erfiaa)) - e^'^' ^ Hn^,{aa) * H ,_,{aa) ^^^^^ 
"n! I 2 (n — ky. ) 



k=0 



2, f 2 

where erf(x) = / is integral of probability fTj. 



VIII Conclusions 

In all problems considered above the attempt to describe the tunneling process of 
particle through barrier with its oscillations between wells in double-well symmetric 
potential was made on the basis of such main parameters as period T of parti- 
cle oscillations between wells, transmission coefficient D through barrier, reflection 
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coefficient R from barrier (for squared potential). The transmission coefficient D 
through barrier is found with consideration of particle, which is incident upon the 
barrier having asymptotic velocity hk/m and is initially determined for continuous 
energy spectrum (see [HI El)- One can obtain all considered above parameters and 
describe the tunneling of particle through barrier from found energy eigenvalues of 
considered systems. If the transmission coefficient is small, then it can be calculated 
in WKB approximation for discrete energy spectrum. Comparison of both methods 
realized for squared potential shows that the values of transmission coefficient, cal- 
culated by these methods, are equal with accuracy to normalization constant which 
is determined by accuracy of used WKB approximation and equal to 1 for ffist order 
terms in h. Therefore, the transmission coefficient of particle through barrier be- 
tween two wells in double-well potential is considered only formally, if it was initially 
determined for continuous spectrum, or with accuracy to normalization constant, if 
it is obtained in WKB approximation. 

The dependencies of transmission coefficient when its value is small on another 
parameters are identical for both methods. These dependencies are obtained, if 
these parameters are the oscillation period between wells and largest common divizor 
determined by expression (1.1) when WKB approximation is applied. 

All parameters considered above and describing the tunneling of particle through 
barrier can be obtained on the basis of found energy levels of considered systems. 
For asymmetric forms of potential they are more difficult to found. 

For asymmetric forms of potentials the splitting of energy levels, which occurs 
because of tunneling of particle through barrier and is studied by theory of WKB 
approximation, gives some features, which are absent when potential is symmetric 
(for example, the possibility of particle reflection from barrier along level, which is 
located higher than height of barrier). In case of potential asymmetry the use of 
boundary conditions must be careful enough. Therefore, in present work a consid- 
erable attention is paid to the problem of finding the solutions of energy levels for 
asymmetric potentials. 

The problem of squared potential was studied early in literature (for example, 
see |H1 m El Hm HI inj)- This problem is considered here as one of the simplest 
cases because it can give visual teaching picture for studding the tunneling particle 
behavior. Relative easiness of finding the solutions in comparison with problems 
with rounded off potentials allows to study deeply the process of tunneling on the 
basis of oscillation period and transmission and reflection coefficients. 

It is significantly difficult to obtain solutions for potential, which has rounded off 
form. Using double- well Morse's potential (as an example), the analysis of solution 
existence is studied and energy eigenvalues equation is obtained, the further solution 
of which can be calculated by use of numerical methods. 

The interesting symmetric double-well potential of form + B"^ /x^ is also ana- 
lyzed because the exact analytical solutions exist for it (and one can obtain energy 
eigenvalues En-, period T of oscillations between wells in the explicit form). This po- 
tential is suitable enough for analysis of particle behavior in wells with enough high 
barrier (with enough large depth). Models of one- dimensional motion of particle in 
such form of potential have been studied in the literature [TT] . 

In the next works it is expected to obtain basic parameters determining behavior 
of particle motion in some forms of asymmetric double-well potentials. 
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Figure 4: The values are exact graphic solutions of system (4.15) and the values 
kn are graphic solutions on this system, where the second equation is changed to its 
linear approximation /2(A;„) 
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Figure 5: Symmetric potential of form x'^ + B^/x^ 
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